Quantum state reconstruction via continuous measurement 
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We present a new procedure for quantum state reconstruction based on weak continuous measure- 
ment of an ensemble average. By applying controlled evolution to the initial state new information 
is continually mapped onto the measured observable. A Bayesian filter is then used to update the 
state-estimate in accordance with the measurement record. This generalizes the standard paradigm 
for quantum tomography based on strong, destructive measurements on separate ensembles. This 
approach to state estimation can be non-destructive and real-time, giving information about observ- 
ables whose evolution cannot be described classically, opening the door to new types of quantum 
feedback control. 
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The control of quantum mechanical systems is finding 
new applications in information jjrocessing tasks such as 
cryptography and computation JJ . Experimental recon- 

I struction of a quantum state is thus essential to verify 
preparation, to detect the presence of errors due to noise 

] and decoherence, and to determine the fidelity of con- 
trol protocols using process tomography. Moreover, real- 

' time "state estimation" may allow improvement of pre- 
cision metrology beyond the standard quantum limit , 

[ with the possibility of active control through closed-loop 
feedback protocols 0]- In addition, measurement of the 

^ quantum state can provide information about unknown 
or nontrivial dynamics, such as those arising in the study 

' of quantum chaos. Laboratory demonstrations of state 
reconstruction are numerous and span a broad range of 

[ physical systems, includinglight fields 0, molecules 0, 
ions joL atoms spins s, 9|, and entangled photon 

[ pairs jldi . 

In this letter we consider a new protocol for quantum 
state reconstruction based on continuous, weak measure- 
ment of a single observable on a single ensemble of iden- 
tically prepared systems. The ensemble is driven so that 
each member undergoes an identical, carefully designed 
dynamical evolution that continually maps new informa- 
tion onto the measured quantity. This is in contrast to 
the standard paradigm for quantum state reconstruction 
based on strong and therefore destructive measurements, 
often of a large set of observables performed on many 
copies of the unknown state. Our weak measurement 
approach has a number of possible advantages in situ- 
ations that lend themselves naturally to working with 
ensembles. Strong measurements on ensembles are in- 
efficient because only a single observable can be mea- 
sured after each preparation and the information gained 
about the observable is extracted independent of the re- 
quired fidelity. By contrast, weak measurements can be 
optimized to obtain just enough information to estimate 
the density matrix to some required fidelity, in real time, 
and with minimal disturbance of each member. One can 



then imagine using the extracted information to perform 
closed-loop feedback control based on knowledge about 
the entire quantum state, rather than on conventional 
"state estimation" of Gaussian random variables that 
evolve according to classical dynamics, as considered to 
date • Our procedure is broadly applicable in systems 
where continuous weak measurement tools have been de- 
veloped, such as nuclear magnetic resonance in molecules 
a and polarization spectroscopy in dilute atomic vapors 
12| , but where noise and decoherence limits the ability to 
perform strong measurements regardless of the amount 
of signal averaging. 

To perform quantum state reconstruction one consid- 
ers a set of measurements {M^'^}, each of which has a 
set of outcomes rn^j'\ {j = 1,2,... , jmax}- The set is said 
to be "informationally complete" if for a given state po 
the set of probabilities {pf^} of the measurement out- 
comes can be inverted to determine po- In our proto- 
col the probabilities are assigned using a single ensemble 
Pn = Po^ J whose dynamical evolution is driven in a 
known fashion and monitored by a probe that measures 
the sum of the identical observables {0} on each mem- 
ber. Due to the central limit theorem the measurement 
record of this probe has the form 

M{t) = N{0), + AM{t), (1) 

where (O)^ is the quantum expectation value at time i, 
and AM{t) is a Gaussian white noise process with vari- 
ance (7^ — l/nAt for measurement strength k and de- 
tector averaging time At. In principle a measurement 
of the collective observable N {O)^ leads to backaction 
on the collective many-body state and can cause indi- 
vidual members of the ensemble to become correlated 
.13j . Such correlations influence the outcome of future 
measurements and greatly complicate the task of recon- 
structing the initial state po- Additionally, the gain from 
performing such quantum limited measurements is small, 
as the majority of the information about the state of indi- 
vidual ensemble members has already been extracted by 
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the probe prior to reaching the quantum hmited regime. 
We thus restrict our considerations to cases where the 
measurement uncertainty, averaged over the total mea- 
surement time T, is large compared to the intrinsic quan- 
tum uncertainty (projection noise) of the collective ob- 
servable, 1/ kT > NAO^, and backaction onto the collec- 
tive state is insignificant. Experimentally this is also the 
most common situation. Of course a sufficient measure- 
ment signal-to-noise ratio must still be available to recon- 
struct the state of an individual member of the ensemble. 
This requires A'^ >> 1 so that the quantum backaction 
associated with information gain is distributed uniformly 
among the entire ensemble, with negligible disturbance of 
any single member state. 

The goal is to invert the measurement history, Eq. Q , 
to determine po- As we wish this procedure to be in- 
dependent of pq, it is most convenient to work in the 
Heisenberg picture and express (O)^ — Tr [0(<)/9o] = 
(0(t)|/9o), where in the second equality we have writ- 
ten the trace as an inner-product between "superopera- 
tors" .14]. We coarse grain over the detector response 
time At, such that {Oi\ = {0{t)\dt/At, obtaining 

a discrete measurement history time-series {Mi}, with 
Mi = N (Oijpo) + (tW, where now the measurement op- 
erators {Oi} are determined in advance by the known 
dynamics, and where is a Gaussian random variable 
with zero mean and unit variance. This equation recasts 
the reconstruction problem as a stochastic linear estima- 
tion problem for the underlying state po- 

In order to reconstruct the state from the measurement 
time-series, the set of measurements operators {Oi} must 
be informationally complete, spanning the space of den- 
sity operators, i.e., the dynamics must map the initial 
measurement to all possible (Hermitian) measurements. 
This is best achieved by introducing an explicit set of con- 
trol parameters, with a time dependent series of Hamil- 
tonians {Hi}. We require that this set generate the Lie 
algebra for SU{d), where d is the dimension of the space; 
the system must be "controllable". Generally, the evo- 
lution will have some associated decoherence that will 
degrade the measurement. Including these terms, the 
evolving measurement operators can then be expressed 
in terms of the base observable as {0{t)\ = (0| St where 

St ~ Texp dt' Lf with Lt the generator of the dy- 
namics, and T the time ordering operator. To simulate 
this evolution we consider a time scale H over which d 
changes negligibly. The semigroup property then allows 
us to approximate iSt+^t = e*'^''^*5t which can be nu- 
merically calculated given a system of reasonable size, 
d < 100. 

A Bayesian filter determines how our knowledge of po 
is updated due to a measurement history {Mi}, 



P{po\{M,}) = AP{{M,}\po)P{po). 



(2) 



distribution and P{po) contains the prior information, in- 
cluding the fact that po is a valid density matrix (ie. has 
trace one and is positive). P{{Mi}\po) is the conditional 
distribution, which contains the information gained dur- 
ing the experiment. The conditional distribution thus 
quantifies how well the measurement performs. Due to 
the Gaussian measurement statistics, this distribution 
has the form, 



P({M,}|po)«exp[- {6p\TZ\6p)]. 



(3) 



The superoperator TZ is the covariance matrix for the 
measurements and Sp = pQ — p is the difference between 
the prepared state and the maximum likelihood estimate 
of this state given the measurements, equivalent to the 
least squares estimator for a Gaussian random variable, 

^ = ^Eio^)(0'i' Ip) = 4e^^^^"'io^)- (4) 



This evolving covariance matrix generalizes the classi- 
cal update rule discussed in 11] . For systems beyond 
spin-1/2, full state reconstruction requires information 
about higher moments of observables, whose evolution 
is fundamentally quantum mechanical. The conditional 
probability distribution has entropy 



5=-log7e=-^logA„ 



(5) 



Here A is the normalization constant for the posterior 



where Xj are the eigenvalues of the covariance matrix, 
corresponding to the inverse of the variances of Eq. (O 
along its primary axes. Qualitatively y/Xj is the signal- 
to-noise-ratio with which we can extract a measure of 
one of the eigen-operators of TZ from the measurement 
record. This entropy thus provides a collective measure 
of the information gained about all parameters, indepen- 
dent of the initial state and any prior information. To 
obtain an accurate reconstruction we need to optimize 
the entropy and any additional costs over the free pa- 
rameters (controls). 

Given the measured information, one estimates the 
quantum state using the mean of the Gaussian condi- 
tional distribution, as given by the least squares fit, Eq. 

Because this does not take into account the prior 
information, the estimated state p is not automatically 
a valid density matrix. To correct this, we must enforce 
Tr [p] = 1 and positivity. The trace condition is ensured 
by adding a pseudo- measurement of the trace Mq = I /d 
which has variance ctq ~ 0. To enforce positivity one 
could solve for the closest positive state to p using convex 
optimization |l5| . Alternatively one can get a reasonable 
and much simpler estimate by setting the negative eigen- 
values of p to zero, and renormalizing to give Ppos, suf- 
ficient when the estimated state has good fidelity. This 
procedure is used in the example below. 

As a concrete demonstration of the power and versa- 
tility of our method we consider the reconstruction of 
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the quantum state associated with the total spin-angular 
momentum of an ensemble of alkali atoms, in our specific 
example the F = 3 or F = 4 hyperfine manifolds of the 
65'i/2 ground state of ^-^^Cs. The number of parameters 
needed for reconstruction are then (2F +1)^ — 1, giving 
48 and 80 components respectively. Consider a cloud of 
cesium atoms prepared in identical states po and cou- 
pled to a common, linearly polarized probe beam tuned 
near the Dl (GS'i/a ^ P1/2) or D2 (651/2 ^ P3/2) res- 
onance [1^ . Information about the atomic spins is ob- 
tained by measuring the Faraday rotation of the probe 
polarization, which provides a continuous measurement 
of the spin component along the direction of propaga- 
tion, O = Fz- Shot noise in the probe polarimeter gives 
rise to the fluctuations W which limit the measurement 
accuracy. 

In the regime of strong backaction onto the collective 
spin state, such measurements have been used to gen- 
erate spin squeezed states l3l [l3l| . and to perform sub- 
shot noise magnetometry [j, lllj . In the regime of neg- 
ligible backaction that is of interest here. Smith et al. 
continuously monitored the Larmor precession of spin in 
an external magnetic field, and observed a series of dy- 
namical collapse and revivals due to a nonlinear term in 
the spin Hamiltonian While this nonlinear collapse 
limits the observation window of a quantum nondemo- 
lition measurement, it also allows for full controllability 
of the atomic spin so that in principle one can recon- 
struct the input quantum state according to the proce- 
dure described above. The nonlinearity results from the 
AC Stark shift caused by off-resonance excitation of the 
Dl or D2 transition, and is directly proportional to the 
excited state hyperfine splitting. Off-resonance excita- 
tion also introduces a small but unavoidable amount of 
decoherence due to photon scattering. Quantum state 
reconstruction requires a large enough nonlinearity to 
generate dynamics that cover the entire operator space 
before decoherence erases information about the initial 
state. This, therefore, favors large excited state hyper- 
fine splittings. 

To control the system we apply a time dependent mag- 
netic field. The overall Hamiltonian, including the non- 
linear AC Stark shift induced by an x-polarized probe, is 

m 

i/(t)-gF//B(B(i) + Bo)-F + /3ft7F2 (6) 

where B(i) is the control field and Bq represents any 
background field that might be present in an experiment. 
In the nonlinear term we have factored out the scatter- 
ing rate 7 for a transition with unit oscillator strength, 
and introduced the ratio /? between the timescales for co- 
herent evolution and decoherence due to optical pump- 
ing. In this example we explore two regimes: a probe 
detuned from the D2 transition by much more than the 
excited state hyperfine splitting {(3 = 0.81), and a probe 
tuned halfway between the two excited hyperfine states 
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FIG. 1: Quantum state reconstruction of the "cat" state de- 
scribed in the text, for the F = 3 hyperfine ground state 
of ^■'■'Cs probed on the Dl transition. The figure shows the 
fidelity of reconstruction as a function of the measurement 
signal-to- noise ratio. Large nonlinearity and little decoher- 
ence allows good fidelity reconstruction in a single run, as 
indicated by the dashed line where a fidelity T = 0.95 is 
reached already at a very modest SNR of 30. The inset shows 
how the reconstruction degrades due to errors in the control 
field, at a constant SNR of 30. For a reasonable 1% control 
field error the fidelity drops to = 0.85. 

(/? = 7.67). Note that the interaction with the mag- 
netic field alone only generates SU{2) rotations which, 
for _F > 1/2, is insufficient to generate the full SU{2F+1) 
algebra. The evolution of the ensemble is governed 
by the master equation Ct[p] = —jr[H{t),p] — 
where p has support only on the ground state of interest 
and all other states have been adiabatically eliminated 
|l7l |. The superoperator includes optical pumping 
within and out of the initial hyperfine manifold. To 
simulate this evolution we construct the £t for some 
choice of scattering rate, background field, and controls 
B(t). The measurement strength k is determined em- 
pirically by the shot-noise limited measurement uncer- 
tainty (7, which we characterize by the signal-to-noise- 
ratio SNR — Mmn^/u, where Mmax = niaxp(F^|p) is 
the maximum signal possible. We account for inhomo- 
geneous values of Bq by averaging over a Gaussian dis- 
tribution corresponding to a standard deviation of 60 Hz 
in the induced Larmor frequency. The duration of the 
simulated measurement is T = 4ms, the coarse graining 
time is 4/is, and the average photon scattering rate is 
7 = lO'^s"^. Finally we simulate the effect of a low pass 
filter by averaging our measurement over a few coarse 
graining time steps. 

The dynamical generation of a complete set of mea- 
surements incurs an unavoidable decoherence cost and so 
must be done as efficiently as possible. To this end we op- 
timize based on a cost function consisting of the entropy 
S{TZ) (Eq. and any additional control costs C(B). In 
this case the additional costs include the degradation in 
reconstruction due to loss of field control at large am- 
plitudes, and inability to rapidly change the fields. We 
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minimize these costs over all possible time dependent 
magnetic fields B(i). This is done sequentially, first re- 
stricting the magnetic fields to optimize the control costs 
C(B), and then optimizing the entropy subject to these 
restrictions. The magnetic fields is restricted to be in the 
X — y plane, with magnitude |B| = 5, such that the Lar- 
mor frequency is lob = HsB/h — 15 kHz. Additionally 
we specify the field at only n = 50 independent times and 
smoothly interpolate between them to ensure slow vari- 
ation. The only free parameters are thus a set of n = 50 
indepedent angles. Optimizing the entropy subject to 
these constraints, we find that the landscape has many 
local minima, precluding the use of purely local search 
techniques. Instead, we use a one dimensional global 
search, where we iteratively optimize one of the n = 50 
independent angles, holding the others fixed. This pro- 
cess is repeated until all of the angles are globally sta- 
tionary, within some tolerance, assuming the others are 
held fixed. This procedure is suboptimal, but converges 
reasonably well. 

The results of simulated reconstructions are shown 
in Figs, n and |21 Given an initial preparation in the 
F ~ 3 hyperfine manifold, in the "cat state" 1-00) = 
(|m 3) + |to = — 3))/-\/2, the fidelity of the reconstruc- 
tion, J-' — {ipol ppos |'0o)i averaged over 1000 noise realiza- 
tions, is plotted versus simulated SNR. Increasing SNR 
clearly results in a better reconstruction. The parame- 
ters needed to attain good fidelity for reconstruction of 
the F = 3 ground state using the Dl transition (Fig.nj 
appear to be well within the reach of current experiments 
at « .95, SNR = 30 Even with a possible 1% un- 
certainty in the control fields, a fidelity of = .85 should 
be possible (Fig.Qlinset). A different regime is illustrated 
by reconstruction of a spin F — A using the D2 resonance. 
This turns out to be infeasible with a single measurement 
run performed on a single ensemble, even assuming very 
large SNR (Fig. EJl. This is because more parameters 
must be reconstructed, 80 versus 48, and because the 
decoherence incurred in conjunction with the nonlinear- 
ity is larger due to the smaller hyperfine splitting of the 
P3/2 manifold. It is still possible, however, to obtain a 
high fidelity reconstruction if we combine the measure- 
ment records from multiple independent runs that each 
start with a fresh ensemble and explore operator space 
in different ways. 

We have presented a new protocol for quantum state 
reconstruction based on continuous measurement of an 
ensemble of N members and demonstrated our proce- 
dure through a simulated reconstruction of a spin J via 
polarization spectroscopy of a gas of cold atoms. The 
reconstruction technique is nondestructive and exploits 
classical estimation theory, providing a starting point for 
consideration of more complex applications of quantum 
control tasks such as quantum feedback. In the future 
work we plan to improve our optimization procedure for 
robustness in control-parameter uncertainty and exam- 
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FIG. 2: Quantum state reconstruction for the F = 4 hyper- 
fine state probed on the D2 transition, where the nonlinearity 
is weaker than for the case in Fig. Q The curves correspond 
to an increasing number of independent runs, which allow in- 
creasingly complete exploration of the operator space. From 
lowest to highest fidelity we use: 1 run (blue) , 2 runs (green) , 
3 runs (red), 4 runs (cyan), 5 runs (magenta), and 6 runs 
(yellow) . 

ine global search procedures such as convex optimization 
(isf . The tools developed here should provide new av- 
enues for real-time state quantum estimation that allow 
us to explore the dynamical generation of nonclassical 
features, such as entanglement. This is of particular in- 
terest for meso scop ic systems whose classical description 
exhibits chaos uS]. 
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